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Abstract 

We consider the equations of state for systems of particles with exclusion statistics 
in a harmonic well. Paradygmatic examples are noninteracting particles obeying ideal 
fractional exclusion statistics placed in (i) a harmonic well on a line, and (ii) a harmonic 
well in the Lowest Landau Level (LLL) of an exterior magnetic field. We show their 
identity with (i) the Calogero model and (ii) anyons in the LLL of an exterior magnetic 
field and in a harmonic well. 

Keywords: exclusion statistics, equation of state, harmonic potential, Calogero model, 
anyons, lowest Landau level 

^Permanent address: Medical Radiological Research Center, Obninsk, Kaluga Region 
249020, Russia. 

^Address from 1 August 1996: Institute of Physics, University of Oslo, P.O. Box 1048 
Blindern, N-0316 Oslo, Norway. 

■^Unite de recherche du CNRS associee aux Universites Paris 11 et Paris 6. 



1 



Haldane's definition of fractional statistics based on a generalized exclusion principle 
IQ when applied locally in phase space results in single state statistical distributions for 
^^^ea/ fractional exclusion statistics (IFES) 0, ^. IFES first appeared in the litterature for 
anyons in the lowest Landau level (LLL) of an exterior magnetic field, i.e. for a gas 
of particles with a degenerate one body spectrum. On the other hand, in one dimension, 
IFES can be modelled with inverse square interactions 0. The relevant systems are the 
Sutherland model (particles on a circle) and the Calogero model (particles on a line 
in a harmonic potential) 0. 

To derive the equation of state and the virial expansion of interacting systems mod- 
elling IFES 1^, 0, one can start from a system placed in a harmonic potential and then 
use a thermodynamical prescription to obtain the equation of state of the original sys- 
tem in the infinite box. The harmonic potential is then referred to as a "long distance 
regulator". This procedure, originally proposed for anyons [10], has been put on general 



ground in JTI], |]. In particular, this is precisely in this way that the equation of state 



for anyons in the LLL was originally derived Q (for several species of particles see [T^.) 
The same procedure can be applied to Id particles with an inverse square interactions: 
virial coefficients for this model have been derived [|| starting from the energy levels of 
the Calogero model in a harmonic well and using an adequate Id thermodynamic limit 
prescription. The latter result was also verified using the equivalence of the Sutherland 



model to a system of free particles obeying IFES |T3|. 

In this note, we address the question of the equation of state for exclusion statistics 
in a harmonic potential, without considering the transition to a system in the infinite 
box limit. We refer to this C3iSG clS db "physical" harmonic potential. What we have in 
mind are mesoscopic quantum dots systems where finite volume effects are supposed to 
be described by a harmonic well. Since there is no well defined volume in this system, the 
equation of state (or the potential Q) depends on the particle number N, the harmonic 
potential frequency u, and other parameters as well. A similar form for the equation 
of state have already been discussed in the Thomas-Fermi approximation for Id IFES 
particles in an external potential [l^ . 
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We start by assuming a certain temperature scaling for the one-particle partition 
function of non interacting IFES particles. In this context, the equation of state is valid 
both for particles in an infinite box and in a physical harmonic well. As examples, we 
discuss non interacting particles in (i) a Id harmonic well and (ii) in a 2d harmonic well 
in the LLL of an external magnetic field. Both these cases correspond to a constant 
density of levels in energy. We compare these equations of state with those for (i) the 
Calogcro model and (ii) anyons in a harmonic well in the LLL of an external magnetic 
field. We discuss more generally a specific form of the equation of state for IFES with a 
constant density of states for systems both in a box and in a harmonic well. Introducing 
the effective volume occupied by the gas in a harmonic well at a given temperature, we 
obtain finally an equation of state in a physical harmonic well quite similar to those in a 
box. 

Ideal fractional exclusion statistics can be defined by the single-state grand partition 
function 



where £j is the energy of the state i and Xi is the Gibbs factor Xi = e^^^~'^^\ (3 — 1/k^T. 

The statistics parameter g = ^ corresponds to Bose and = 1 to Fermi statistics. The 
distribution function n{xi) = Xi{d/dxi)hi^{xi) is connected with ^{xi) by the bilinear 
relation 




(1) 



1 



1 



(2) 



n{xi) 



+ 9- 



Expanding in powers of Xi yields 



k=l 

By summing ln^(a;j) and n{xi) over i, 




(3) 



QO CO 



(4) 



fc=l k=l 



{z = & 



is the fugacity), with the cluster coefficients 



(5) 
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ZliP) = Z^i is the one-particle partition function. 

Without specifying in detail the system at this stage, we assume that the partition 
function Zi{(3) scales with the inverse temperature (3 as 

This scaling is relevant for systems with a gap £0 in the single particle energy spectrum, 
with one-particle partition function Zi{(3) = e~^^" Z[{i3). This factorization will indeed 
materialize in the thermodynamic limit for various physical systems as we will see below. 
Taking into account (H), the cluster coefficients (^) become 

h = ^e-'^^^^Zm . (7) 

One deduces from (^) the "virial expansion" 

k=l 



with the dimensionless "virial coefficients" 

Ai = l, A, = ^{2g-1), A, = [{A-'-'-2.3-'-')+g{g-l){A-'-3-%... (9) 

We now specialize on the particular case 5 = 0, which can be completely analyzed. 
This case is of particular interest since it implies a scaling for Zi{j3), which means a 
constant density of states in energy. One can already see from (H) that 5 = is special 
since it implies that ^3 does not depend on the statistical parameter g. 

Comparing the second expansion in (Q), with the cluster coefficients (J^), with (||), one 
can write 

N = Z[{l3)\naz'), (10) 

where z' = ze~^^° . Regarding Eq. (p!OD as determining z' as a function of N , we obtain 
from (H) 

_djn^idzv^ 

dN z ON ' ^ ' 
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To calculate dz\N)/dN (or dN{z')/dz'), we use (|T0|) and ([I])-®. Noting that 

{l/N{z'))dN{z')/dz' = n{z')/z', we find an equation of state 

_d(3n_ N / 1 \ 
" 'W ~ Z[{P) [e^/K(P) - 1 + ^) ' 

and upon integration, finally, 

~^ Jo Zm (e^'/^i('3) _ 1) + 2^^7(;^ ■ (13) 
Expanding this, we obtain the virial expansion 

where are the Bernoulli numbers [B^ — — vanishing for k odd. It appears 

that the statistical parameter g only enters the equation of state via the second virial 
coefficient 



Let us now examine how the equations of state (H) and (|T3D can be relevant to both 
systems in a box and in a physical harmonic well. 

We first consider a gas of free (spinless) particles, with generic dispersion law = 
+ O'V'^ occupying a box of volume in d dimensions. In the thermodynamic limit the 
one-particle partition function is 

ZJl,Q\ - r-^0e, T{l + d/a)V 

^'^"'^^ ~ ' (2v/?)'^r(l + ld){akf3)-'^/'^ ' ^^^^ 
satisfying (|]), with 6 = d/a — 1. In this case the virial expansion @ becomes the usual 
virial expansion for a system in a box with pressure P given by 

oo 

pP = J2akp', P = N/V, (16) 

k=l 



with the (dimensional) virial coefficients 



where are given in (^). Note that the virial coefficients (0) coincide with those 
obtained in |T^ for £0 = showing that the presence of a gap in the particle dispersion 



does not affect the equation of state. 



The constant density of states case, 5 = d/ a — 1 = describes for example chiral 
particles on a line with linear dispersion with a gap, 

e = eo + vp , p > . (18) 

The equation of state is (|T3D (or (0)) with 

zm = ^ • (19) 

27ivp 

Let us now turn to systems in a "physical" harmonic well. We restrict to cases with a 
constant density of states, i.e. we consider noninteracting (nonrelativistic) particles that 
occupy single particle levels 

(i) in a Id harmonic well V{x) = ^moj'^x^ ; 

(ii) in the LLL of an external magnetic field i? in a 2d harmonic well. 

In a quantum dot language, where uj is small (for example with respect to the cyclotron 
frequency), but non vanishing, the harmonic well is supposed to encode the finite size 
effects of the sample. 

The one particle energy levels are 

e^ = eQ + ^w , ^ > , (20) 

where for the Id harmonic model 

1 , , 

eQ = -UJ , w = uj , (21) 

whereas for the 2d LLL harmonic model, one has 



eo = \/t^c+^^5 = + uj'^ - ujc , (22) 
or, to leading order in uP' I<jj\ 



,2 



£0 = tUc , w =^ . (23) 
It follows from (^OD that the one-particle partition function is 

-A;/3eo p-fc/3eo 

zm) = = + ^^^^^) + ^[(^/^-)']) • (24) 



1 - e-^^^ kl3vj 
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Thus, to leading order in k(3tu, the partition function (p^ scales as in (|), with 5 = and 

Z[{P) = 1/Pzu . (25) 

This leads to the equation of state (^),(0). If one considers the correction terms 0{k(3w) 
in (p4D, they might lead to corrections to the virial coefficients of very high order, and 
above. However, if the virial expansion converges, these corrections are negligible inside 
the radius of convergence. 

It is interesting to note that for Id IFES particles in a harmonic well, the equation 
of state ([T3| ) ,([T^ ) with (^) and (pT]) coincides with that obtained in the Thomas- Fermi 
approximation . We stress that the present derivation only uses the condition /3a; <^ 1 
when the discreteness of energy levels becomes inessential. 

Interacting systems modelling IFES are known to be the Id Calogero model and 2d 
anyons in the LLL of an external magnetic field. An harmonic well is added to these 
systems to define a proper thermodynamic limit. It would be interesting to investigate 
what happens if the harmonic well becomes now physical. The iV-body energy spectrum 
is given by [|] 

Y,n,e, + \N{N -l)wg , (26) 

where ni are nonnegative integers with the constraint Yl,i = N, and where ee is given 
by (^). For anyons, (7 = a G [0, 1] is the anyonic statistical parameter; for the Calogero 
model, g = a >0 defines the singular two-body potential a{a — l)/{xi — Xj)'^ and specifies 
the coinciding point (short distance) behavior of the A^-body wave function, \E' 

9j S I OC 2 3y j I ^ • 



The A^-body partition function corresponding to (|26D is 



;„ = e-4^-<"-)- n ■ (27) 

n=l ^ 



Remarkably, the cluster coefficients obtained form ( P7[) are 



hu = e-'P^-^{l + 0{Pw)), (28) 
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having, to leading order in Pzu, the form (|^) with 5 = and Z[{P) given by (|25D. As above, 
the terms of order 0{(3w) in (pSf ) are neghgible inside the radius of convergence of the 
virial expansion. It follows that both equations of state of the Calogero and LLL-anyon 
model in a 2d harmonic well are given by ([T^),(|T^) with Z[{j3) determined by (|TB|), i.e. 
they are identical to the equation of states for noninteracting particles in a Id harmonic 
well or in a 2d LLL harmonic well. 

This happens to be general feature of equations of states in a physical harmonic well 
for various physical systems modelling IFES with a constant density of levels: only the 
second virial coefficient depends on the statistics parameter, in agreement with the general 
statements for free IFES particles in a box | |13| ]. 

As a specific example, consider the equation of state for a model of a chiral field on 
a circle proposed in [jl6|. This model was constructed by mapping the second quantized 
LLL anyon model on a circle. The energy levels of this model are given by ( PB| ) with the 
identification 

2uJc V 

where V is the length of the circle. If the harmonic potential is assumed to materialize in 
the anyon droplet of radius V/27T by the action of an electric field, the velocity v can then 
be interpreted as the drift velocity E/B on the edge (the velocity of the edge excitations), 
where the electric field on the edge is E = {m/ e)uj'^R. The thermodynamic limit is 
understood as ^ 0, ^ oo whereas v is kept fixed. The spectrum (^), altogether 
with the identification (pQ]), yield the cluster coefficients (|^) with 5 = 0, eo = {.Z[{[3) 
is given by (0)). It follows that the equation of state is (PTB]) , (]T^) , that is those of free 
chiral Id IFES particles with dispersion (|18]). This conclusion is consistent with those of 
T7| , namely the model of a chiral field on a circle fl^ admits an interpretation in terms 



of IFES. 



We finally note that (|T3D (or (14)) for a physical harmonic well, with Z[[[3) given 
by (^), can be viewed as an equation of state connecting the average pressure to the 
average particle density in the harmonic well. The local pressure -P(x) in a slowly varying 
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(i- dimensional harmonic potential is defined as 

= J (i'^X;3P(x) . (30) 

The average pressure is then 

(3{P)^(^^^ Jd'^pPi^), (31) 

where Ay = \j2'n(3/m is the thermal de Broglie wavelength. The local particle density 
p(x) being normalized as / d^^p{:s) = N, one can define the average density as 

{p) ^ 1^^^' I d'M^) . (32) 
Introducing the effective volume occupied by the gas at temperature T by 

V..-(^)\ (33) 

we have Kfr ~ -Rfg, where Res is determined by 

Imu'R'^^ ^ T . (34) 
It follows that (p!4D can then be rewritten as a virial expansion in a box of volume Ves (cf 

(0), my- 

l + -(2^-l)/3^Kff(p) + E^^^(/5^Kff)'=(p)'j • (35) 

For noninteracting particles in a Id harmonic well and for the Calogero model, j3toVcs = 
Xt, whereas for noninteracting LLL particles in a 2d harmonic well and for LLL anyons 
in a 2d harmonic well, jSwVes = 27i/eB. For the Calogero model, the equation of state 
(p5| ) was conjectured in Note that for the LLL anyon model, the condition 1 <^ 2(3uJc 
guarantees that the particles, when confined in the quantum dot of size i?efr given by (^l)), 
do not reach the second Landau level. 
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